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Background: Studies of fission dynamics, based on nuclear energy density functionals, have shown that the 
coupling between shape and pairing degrees of freedom has a pronounced effect on the nonperturbative collective 
inertia and, therefore, on dynamic (least-action) spontaneous fission paths and half-lives. 


Purpose: To analyze effects of particle-number fluctuation degree of freedom on symmetric and asymmetric 
spontaneous fission (SF) dynamics, and compare with results of recent studies based on the self-consistent Hartree- 
Fock-Bogoliubov (HFB) method. 


Methods: Collective potentials and nonperturbative cranking collective inertia tensors are calculated using the 
multidimensionally-constrained relativistic mean-field (MDC-RMF) model. Pairing correlations are treated in 
the BCS approximation using a separable pairing force of finite range. Pairing fluctuations are included as 
a collective variable using a constraint on particle-number dispersion. Fission paths are determined with the 
dynamic programming method by minimizing the action in multidimensional collective spaces. 


Results: The dynamics of spontaneous fission of 7°*Fm and ?°°Fm are explored. Fission paths, action integrals 
and corresponding half-lives computed in the three-dimensional collective space of shape and pairing coordinates, 
using the relativistic functional DD-PC1 and a separable pairing force of finite range, are compared with results 
obtained without pairing fluctuations. Results for °°Fm are also discussed in relation with those recently obtained 
using the HFB model. 


Conclusions: The inclusion of pairing correlations in the space of collective coordinates favors axially symmetric 
shapes along the dynamic path of the fissioning system, amplifies pairing as the path traverses the fission barriers, 


significantly reduces the action integral and shortens the corresponding SF half-life. 


PACS numbers: 21.60.Jz, 24.75.+i, 25.85.Ca, 27.90.+b 


I. INTRODUCTION 


Important advances have recently been reported in mi- 
croscopic modeling of the dynamics of spontaneous and 
induced fission, based on nuclear density functional the- 
ory (DFT) [1]. Within this framework, spontaneous fis- 
sion (SF), in particular, is described by quantum tunnel- 
ing through potential barrier(s) in a multidimensional 
space of coordinates that parametrize large-amplitude 
collective motion. Most calculations of SF lifetimes 
are based on the semiclassical Wentzel-Kramers-Brillouin 
(WKB) approximation for the one-dimensional barrier 
tunneling. The dynamics of the SF process is governed 
by the potential energy surface (PES) as function of the 
collective coordinates, and by the collective inertia along 
the fission path. The path along which the nucleus tun- 
nels is determined by minimizing the fission action inte- 
gral in the multidimensional collective space [2, 3]. 

The PES’s can be computed using the macroscopic- 
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microscopic (MM) model, or a number of self-consistent 
mean-field (SCMF) approaches based on microscopic ef- 
fective interaction or energy density functionals. In 
most recent studies the multidimensional collective iner- 
tia tensor is usually determined using the adiabatic time- 
dependent Hartree-Fock-Bogoliubov (ATDHFB) method 
with the perturbative cranking approximation (neglect- 
ing the contribution from time-odd mean fields and treats 
perturbatively the derivatives of single-nucleon and pair- 
ing densities with respect to collective coordinates), or 
the nonperturbative cranking approximation (the deriva- 
tives with respect to collective coordinates are computed 
explicitly). Although perturbative cranking ATDHFB 
collective masses has extensively been used in SF fission 
half-life calculations [2, 4-7], a number of recent studies 
[8-10] have indicated the essential role of the nonpertur- 
bative cranking ATDHFB approximation to the collec- 
tive inertia for a quantitative dynamic description of SF. 

In the first approximation the effective collective in- 
ertia M x A~?, and the collective potential V œ 
(A — Ao)’, where A is the pairing gap and Ao corre- 
sponds to its self-consistent stationary value. When the 
gap parameter is treated as a dynamical variable, an en- 
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hancement of pairing correlations reduces the effective 
inertia and thus minimizes the action integral S along 
the fission path [11]. A number of studies of SF have 
shown that the coupling of pairing fluctuations with the 
fission mode can significantly reduce the estimated fission 
lifetimes [12-19]. 

In recent studies dynamic fission paths determined 
with the least-action principle have been investi- 
gated using the Hartree-Fock-Bogogliubov (HFB) frame- 
work based on the Barcelona-Catania-Paris-Madrid [20], 
Gogny D1M [20], and Skyrme SkM* [21, 22] energy den- 
sity functionals. The pairing gap parameter has been in- 
cluded as a dynamical variable in the collective space. As 
a result, an enhancement of pairing correlations along fis- 
sion paths and the speedup of SF have been predicted. It 
has also been noted that pairing fluctuations can restore 
axial symmetry in the fissioning system [21, 22], although 
the triaxial quadrupole degree of freedom is known to 
play an important role around the inner and even outer 
barriers both along the static fission path for actinide 
nuclei (Ref. [23] and references therein), and in the dy- 
namic case when the influence of pairing fluctuations is 
not taken into account [9, 10]. 


In Ref. [10] we have used the multidimensionally- 
constrained relativistic Hartree-Bogoliubov (MDC-RHB) 
to analyse effects of triaxial and octupole deformations, 
as well as approximations to the collective inertia, on 
the symmetric and asymmetric spontaneous fission dy- 
namics. Based on the framework of relativistic energy 
density functionals, and using as examples 7°4Fm and 
250Fm, our analysis has shown that the action integrals 
and, consequently, the half-lives crucially depend on the 
approximation used to calculate the effective collective 
inertia along the fission path. While the perturbative 
cranking approach underestimates the effects of struc- 
tural changes at the level crossings, the nonperturba- 
tive collective mass is characterized by the occurrence 
of sharp peaks on the surface of collective coordinates, 
that can be related to single-particle level crossings near 
the Fermi surface, and this enhances the effective inertia. 
In this work we continue to explore the dynamics of SF of 
264Fm and 75°Fm but, in addition to shape deformation 
degrees of freedom, pairing correlations are included in 
the space of collective coordinates. The dynamic (least- 
action) fission paths are determined in three-dimensional 
(3D) collective spaces, and the corresponding SF half- 
lives are computed. Since calculations in the 3D collec- 
tive space with the MDC-RHB model are computation- 
ally very demanding, here we employ the MDC-RMF 
model in which the pairing correlations are treated in 
the BCS approximation. The collective inertia tensor is 
calculated using the self-consistent relativistic mean-field 
(RMF) solutions and applying the ATDHFB expressions 
in the nonperturbative cranking approximation. The ar- 
ticle is organized as follows: the method for calculating 
dynamic fission paths is described in Sec. II; numerical 
details of the calculation, results for the deformation en- 
ergy landscapes, collective inertias, minimum-action fis- 


sion paths and the corresponding half-lives are discussed 
in Sec. III; and Sec. IV contains a short summary of the 
main results. 


II. METHOD FOR CALCULATING DYNAMIC 
FISSION PATHS 


RMEF-based models present a particular implementa- 
tion of the relativistic nuclear energy density functional 
(EDF) framework, which has become a standard method 
for studies of the structure of medium-heavy and heavy 
nuclei [24-29]. As in our previous study of spontaneous 
fission [10], here we employ the point-coupling relativis- 
tic EDF DD-PC1 [30]. Starting from microscopic nu- 
cleon self-energies in nuclear matter, and empirical global 
properties of the nuclear matter equation of state, the 
coupling parameters of DD-PC1 were fine-tuned to the 
experimental masses of a set of 64 deformed nuclei in the 
mass regions A ~ 150 — 180 and A ~ 230 — 250. The 
functional has been further tested in a number of mean- 
field and beyond-mean-field calculations in different mass 
regions. 

For a quantitative description of open-shell nuclei it 
is necessary to consider also pairing correlations. In the 
MDC-RMF model pairing is taken into account in the 
BCS approximation and here, as in Ref. [10], we use a 
separable pairing force of finite range: 


V (r1, r2,r1,r5) = Go (R — R’)P(r)P(r’)=(1- P’), 


(1) 
where R = (rı +r2)/2 and r = rı — r2 denote the center- 


of-mass and the relative coordinates, respectively, and 
P(r) reads 


NI = 


1 —r? /4a? 
P(r) = amar ue (2) 
The two parameters Go = —738 MeV fm~? and a = 
0.644 fm [31] have been adjusted to reproduce the density 
dependence of the pairing gap in nuclear matter at the 
Fermi surface calculated with the D1S parameterization 
of the Gogny force [32]. 

The energy landscape is obtained in a self-consistent 
mean-field calculation with constraints on mass multipole 
moments Q),, = T*Y),, and the particle-number disper- 
sion operator AÑ? = Ñ? — (Ñ)? [33]. In the present 
analysis the Routhian is therefore defined as 


1 z 
E' — ERMF + 5 gO Qan + A` AN? 3 (3) 
Ap 


where Epmr denotes the total RMF energy including 
static BCS pairing correlations. The amount of dynamic 
pairing correlations can be controlled by the Lagrange 
multipliers A2; (T = n, p), [21, 34, 35]. As it has recently 
been shown in a similar study of Ref. [21], the isovector 
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pairing degree of freedom appears to play a far less impor- 
tant role in spontaneous fission as compared to isoscalar 
dynamic pairing. Therefore, the computational task can 
be greatly reduced by considering only dynamic pairing 
with Ag = Azp = A2 as a collective coordinate. 

The nuclear shape is parameterized by the deformation 
parameters 


At 


Pre = SARK 


Qan). (4) 
The shape is assumed to be invariant under the exchange 
of the x and y axes and all deformations S), with even 
u can be included simultaneously. The deformed RMF 
equations are solved by an expansion in the axially de- 
formed harmonic oscillator (ADHO) basis [36]. In the 
present study of transactinide nuclei calculations have 
been performed in an ADHO basis truncated to Ny = 16 
oscillator shells. For details of the MDC-RMF model we 
refer the reader to Ref. [23]. 

The action integral along the one-dimensional fission 
path L is calculated using the expression: 


m ~\/ Mens) Vat) mld, 18) 


S(L) = 


Sin 


where Meg(s) and Veg(s) are the effective collective in- 
ertia and potential along the path L(s), respectively. Eo 
is the collective ground state energy, and the integra- 
tion limits correspond to the classical inner (sin) and 
outer (Sout) turning points defined by: Veg(s) = Eo. 
The fission path L(s) is determined in the multidimen- 
sional collective space by minimizing the action integral 
of Eq. (5) [2, 3]. The spontaneous fission half-life is cal- 
culated as T,/2 = In2/(nP), where n is the number of 
assaults on the fission barrier per unit time [9, 21, 37, 38], 
and P is the barrier penetration probability in the WKB 
approximation 


1 


p= 1 + exp[2S(L)| (6) 

The action integral Eq. (5) and, therefore, the fission 
half-life is essentially determined by the effective collec- 
tive inertia and potential. The effective inertia is defined 
in terms of the multidimensional collective inertia tensor 
M [2, 9, 10, 21, 37, 38] 


Mes(s) = 5° M 


a 


ij a =, (7) 
s ds 
where q;(s) denotes the collective variable as function of 
the path’s length. 
In the present study the inertia tensor is computed us- 
ing the ATDHFB method in the nonperturbative crank- 
ing approximation [8] 


h2 
MS 2 


2095 op 


Ea + Es 


(8) 


where 
Fi Op On Op* On* 
— =UiI—v* + Ut U* — vi e_U* — Viv". 
qi Ogi Ogi Ogi Ogi 


(9) 
U and V are the self-consistent Bogoliubov matrices, and 
p and « are the corresponding particle and pairing density 
matrices, respectively. The derivatives of the densities 
are calculated using the Lagrange three-point formula for 
unequally spaced points [39, 40]. 

The collective potential Vag is obtained by subtract- 
ing the vibrational zero-point energy (ZPE) from the to- 
tal RMF constrained energy surface [9, 10, 21, 41, 42]. 
The fission path is determined in a multidimensional 
collective space using both the dynamic programming 
(DPM) [38] and Ritz [37] (RM) methods . For both 
methods we have considered several possible values for 
the turning points sj, and Sout to verify that the min- 
imum action path is chosen. Since both methods give 
virtually identical results, only those obtained using the 
DPM are included in the presentation. 


III. SPONTANEOUS FISSION OF 7°*Fm AND 
250Fm: PAIRING-INDUCED SPEEDUP 


To study the effect of dynamic pairing correlations 
along fission paths, as in our previous study of SF of 
Ref. [10], we will analyze two illustrative examples: the 
symmetric spontaneous fission of 7°4Fm and the asym- 
metric SF of 7°°Fm. In addition to shape variables, 
here pairing correlations are also considered as collec- 
tive coordinates in the study of fission dynamics. Be- 
cause of computational restrictions and to simplify the 
interpretation of results, the present analysis is restricted 
to a three-dimensional (3D) collective space, defined 
by either (820, G22,A2) (quadrupole triaxial shapes) or 
(820, 830, 2) (quadrupole and octupole axial shapes), 
where the coordinate A2 represents dynamic pairing fluc- 
tuations. The relativistic energy density functional DD- 
PC1 [30] is employed in self-consistent RMF calculations 
of constrained energy surfaces, collective inertia tensors 
and fission action integrals. The height of fission bar- 
riers is sensitive to the strength of the pairing interac- 
tion [43] and, therefore, a particular choice of the pair- 
ing strength may have a considerable effect on fission 
dynamics. As explained above and in Ref. [10], the pa- 
rameters of the finite range separable pairing force were 
originally adjusted to reproduce the pairing gap at the 
Fermi surface in symmetric nuclear matter as calculated 
with the Gogny force D1S. A number of mean-field stud- 
ies based on the relativistic Hartree-Bogoliubov (RHB) 
model have shown that the pairing strength needs to 
be fine-tuned in some cases, especially for heavy nu- 
clei [44, 45]. Since in the present study pairing corre- 
lations are treated in the BCS approximation, we have 
adjusted the strength parameters to reproduce the avail- 
able empirical pairing gaps in Fm isotopes. The resulting 
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FIG. 1. (Color online) Effective collective potential Veg of 


264Fm in the (620,822) plane for A2 = 0 (a), and in the 
(820, A2) plane for 822 = 0 (b). In each panel energies are 
normalized with respect to the corresponding value at the 
equilibrium minimum, and contours join points on the surface 
with the same energy (in MeV). The energy surfaces are cal- 
culated with the relativistic density functionals DD-PC1 [30], 
and the pairing interaction Eq. (1). 


values with respect to the original pairing strength ad- 
justed in nuclear matter (Gp = —738 MeV fm~) are 
G,/Go = 1.21, Gp/Go = 1.14. As in Refs. [9, 21] and 
our previous work Ref. [10], we choose Ep = 1 MeV in 
Eq. (5) for the value of the collective ground state energy. 
This value enables a direct comparison of our results with 
those reported in previous studies, especially in Ref. [21]. 
For the vibrational frequency wo = 1 MeV the num- 
ber of assaults on the fission barrier per unit is 1020-35 
s~! [46]. 


A. Symmetric fission of 7°4Fm 


The first example in our analysis of the influence of dy- 
namical fluctuations in shape and pairing degrees of free- 
dom on fission paths is the nucleus 2°4Fm, for which the- 
oretical studies [47, 48] predict a symmetric spontaneous 
fission decay. The shape degrees of freedom in this case 
are elongation and triaxiality and, therefore, calculations 
of the energy landscape, inertia tensor, and fission paths 
are restricted to the 3D collective space (20, B22, A2). 
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FIG. 2. (Color online) Cubic root determinants of the 


nonperturbative-cranking inertia tensor MF]! (in 10 x 
h?MeV—') of 764Fm in the (20,822) plane for A2 = 0 (a), 
and in the (820, A2) plane for 822 = 0 (b). 


In Fig. 1 we plot the collective potential energy (the 
vibrational ZPE is subtracted from the constrained self- 
consistent mean-field energy) of 7°*Fm in the (G20, B22) 
plane for Az = 0 (a), and in the (29, Az) plane for B22 = 0 
(b). Therefore, panel (a) displays the results obtained 
without including dynamical pairing correlations, and 
the potential can be directly compared to the one shown 
in Fig. 6 (b) of Ref. [10], where we used the relativistic 
Hartree-Bogoliubov (RHB) model to compute the energy 
surface. We notice that the energy landscapes obtained 
with RMF and RHB models are almost identical, and this 
validates the treatment of pairing in the BCS approxi- 
mation in the present analysis. The functional DD-PC1 
predicts an axially symmetric equilibrium state at mod- 
erate deformation (G29 ~ 0.2), and the axially symmetric 
barrier at G29 œ% 0.6 is bypassed through the triaxial re- 
gion, thus lowering the height of the barrier by + 2.5 — 3 
MeV. In panel (b) of Fig. 1 we project the potential en- 
ergy calculated in the 3D collective space that includes 
dynamic pairing, on the (G29,A2) plane. For b22 = 0 
(axially symmetric shape), the potential energy increases 
monotonically with A2 (stronger pairing) at each defor- 
mation 629. The topography of the collective potential 
in the (829, 2) plane is relatively simple. 

The effective inertia that determines the fission action 
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integral is defined in terms of the multidimensional col- 
lective inertia tensor M Eq. (7). The important effects 
related to the exact treatment of derivatives of single- 
particle and pairing densities in the ATDHFB expres- 
sions for the mass parameters were recently analyzed in 
Refs. [9] and [10]. For the three-dimensional space of 
collective coordinates, six independent components de- 
termine the inertia tensor. The inertia tensor can be visu- 
alized by plotting the cubic root determinant |M|!/3, in- 
variant with respect to rotations in the three-dimensional 
collective space [21]. 

In Fig. 2 we plot |M°|!/3 obtained in the nonpertur- 
bative cranking approximation, in the (820, G22) plane for 
A2 = 0 (a), and in the (820, A2) plane for 622 = 0 (b). 
These results correspond to the self-consistent solutions 
for the potential energy surfaces shown in Fig. 1. The 
nonperturbative |M°|!/3, calculated without dynamic 
pairing correlations, displays a rather complex structure 
in the (620, 822) plane (Az = 0), as shown in panel (a). 
In particular, very large values of |M|!/8 are calculated 
in the region of the axial fission barrier. As discussed in 
Refs. [8-10], this is related to single-particle level cross- 
ings near the Fermi surface. The abrupt changes of oc- 
cupied single-particle configurations lead to strong vari- 
ations in the derivatives of densities in Eq. (9) and, con- 
sequently, sharp peaks develop. When dynamic pairing 
correlations are included in the collective space (panel 
(b)), a simple dependence of the nonperurbative |M€|!/3 
on Az is obtained at each deformation 329, consistent with 
the expected relation M œ A~?. We note that the re- 
sults for 2°4Fm, shown in Figs. 1 and 2, are very similar to 
those obtained using the nonrelativistic HFB framework 
based on the Skyrme energy density functional SkM* and 
a density-dependent pairing interaction (cf. Fig. 2 of 
Ref. [21]). 

The coupling between shape and pairing degrees of 
freedom has a pronounced effect on the predicted fis- 
sion paths. As the effective potential increases from the 
self-consistent values because of the enhancement of pair- 
ing, the effective collective inertia is reduced x A~?. 
These two effects determine the minimum-action path 
in Eq. (5). The projections of the 3D spontaneous fis- 
sion path of 7°4Fm on the (629, G22) plane and on the 
(G20, 2) plane are shown in Fig. 3 (a) and (b), respec- 
tively (solid curves). The two-dimensional (2D) path cal- 
culated without pairing fluctuations (Az = 0) is also in- 
cluded for comparison (dash-dot curve). It is very in- 
teresting to note that, while the 2D dynamic path de- 
tours the axial barrier through the triaxial region, the 
extension of the collective space by the pairing degree of 
freedom fully restores the axial symmetry of the fission- 
ing system. The evolution of the pairing strength along 
the axially symmetric fission path is shown in the lower 
panel of Fig. 3. One notices how, in order to reduce 
the collective inertia, the fissioning nucleus favors an in- 
crease in pairing over the static self-consistent solution, 
at the expense of a larger potential energy. Because of 
pairing fluctuations, the corresponding fission action in- 
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FIG. 3. (Color online) Projections of the three-dimensional 
(3D) dynamic path (solid curves) for the spontaneous fission 
of 7°4Fm on the (620, 622) plane for àz = 0 (a), and the 
(820, A2) plane for B22 = 0 (b), calculated using the dynamic 
programming method. The dash-dot-dot curve denotes the 
two-dimensional (2D) path computed without the inclusion 
of dynamic pairing correlations. 


tegral is reduced by about 5 units with respect to the 2D 
path and, consequently, the predicted half-life is almost 
five orders of magnitude shorter in comparison to the 
2D case without the dynamic pairing degree of freedom 
(cf. Table I). This result can directly be compared to the 
one obtained using the the Skyrme energy density func- 
tional SkM* and a density-dependent pairing interaction 
(cf. Fig. 3 of Ref. [21]). In the latter case triaxiality 
is reduced along the 3D fission path because of dynamic 
pairing fluctuations, but the full axial symmetry is not 
restored. This is probably due to the fact that in the 2D 
calculation with the Skyrme functional the triaxial coor- 
dinate reduces the fission barrier height by more than 4 
MeV (less than 3 MeV in the present calculation with 
DD-PC1). A combination of a higher axially symmetric 
fission barrier and/or possibly weaker pairing, prevents 
the full restoration of axial symmetry along the 3D fis- 
sion path of 764Fm. In the case of 74°Pu, on the other 
hand, for which the Skyrme functional SkM* predicts an 
energy gain on the first barrier resulting from triaxiality 
of only 2 MeV, the inclusion of pairing fluctuations leads 
to a full restoration of axial symmetry along the 3D fis- 
sion path between the equilibrium ground state and the 
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FIG. 4. (Color online) Same as described in the caption to 
Fig. 1 but for the nucleus 7°°Fm. 


superdeformed fission isomer (cf. Fig. 5 of Ref. [21]). 


B. Asymmetric fission of 7°°Fm 


In the second example we explore the interplay be- 
tween reflection-asymmetric shapes and pairing degrees 
of freedom, and analyze the asymmetric spontaneous fis- 
sion of 7°°Fm [48]. Since the triaxial degree of freedom 
is particularly important around the inner fission bar- 
rier, and the complete calculation in the four-dimensional 
collective space (820, 822, 630, A2) is computationally 
too demanding, we first analyze the path that connects 
the mean-field equilibrium (ground) state and the iso- 
meric fission state calculated in the (29, 822, A2) collec- 
tive space. The collective potential energy surface of 
250Fm in the (820, G22) plane for Ag = 0, and in the 
(820, A2) plane for G22 = 0, is plotted in upper and lower 
panels of Fig. 4, respectively. The inclusion of the tri- 
axial degree of freedom reduces the inner fission barrier 
height by ~ 2 MeV, and this effect is similar in mag- 
nitude to the case of ?°4Fm considered in the previous 
section. The lower panel displays the projection of the 
potential energy calculated in the 3D collective space on 
the (29,2) plane and we notice that, for 822 = 0, the 
energy increases monotonically with Ag at each value of 
the axial deformation parameter (29, with a pronounced 
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FIG. 5. (Color online) Same as described in the caption to 
Fig. 2 but for the nucleus 7°°Fm. 


fission barrier around B29 ~ 0.55. 

The deformation dependence of the nonperturbative 
collective inertia tensor is displayed in Fig. 5, where 
we plot the cubic root determinants |M°|!/3 in the 
(820, G22) and (629, A2) planes. The global deformation 
dependence of |M©|1/% is similar to the one calculated 
for 264Fm and shown in Fig. 2, that is, |MC|!/> dis- 
plays strong variations in the (6820, 622) plane for A2 = 0, 
and pronounced peaks generated by single-particle level 
crossings near the Fermi surface appear in the region of 
the fission barrier. By including the dynamic pairing de- 
gree of freedom, one finds that |M|!/? decreases as 2 
increases at each deformation (29. 


The projections of the 3D dynamic path determined in 
the (620, G22, A2) collective space is shown in Fig. 6 (solid 
curves). The minimum-action path connects the inner 
turning point and the isomer minimum at B99 ~ 0.95. 
The 2D path calculated in (620, 822) collective space is 
also included (dash-dot red curve) for comparison. Even 
though the 2D dynamic path does not extend very far 
in the triaxial region, the triaxial shape degree of free- 
dom is important in the calculation of the fission action 
integral, similar to the result we obtained with the RHB 
model in our previous study of fission dynamics (cf. Fig. 
12 of [10]). However, since triaxiality gains only ~ 2 
MeV in energy on the first barrier, the inclusion of dy- 
namic pairing fluctuations fully restores axial symmetry 
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FIG. 6. (Color online) Projections of the three-dimensional 
(3D) dynamic fission path (solid curves) of 7°°Fm on the 
(820, 822) plane for A2 = 0 (a), and the ((20,A2) plane for 
B22 = 0 (b), The dash-dot (red) curve denotes the two- 
dimensional (2D) path computed without the inclusion of dy- 
namic pairing correlations. The minimum-action paths con- 
nect the inner turning point and the fission isomer. 


in the fissioning system (upper panel of Fig. 6). Pair- 
ing is enhanced with respect to the static solution along 
the axially symmetric path and, consequently, the action 
integral in the interval between the inner turning point 
and the isomeric state decreases from 23.06 (2D dynamic 
path) to 14.90 along the 3D fission path. 

Since the triaxial shape degree of freedom is sup- 
pressed, that is, it does not contribute to the action 
integral in the dynamic case when pairing fluctuations 
are included, we can analyze the SF decay of 7°°Fm 
in the restricted 3D collective space with coordinates 
(820, 830, 2). In Fig. 7 we display the collective po- 
tential of 7°°Fm in the (20, 830) plane for Az = 0 (a), 
and in the (820, A2) plane for 639 = 0 (b). The mean- 
field equilibrium (ground) state is predicted at moder- 
ate quadrupole deformation 629 + 0.3, and the isomeric 
minimum at B29 % 0.95. The nucleus remains reflection 
symmetric through the entire region of quadrupole de- 
formations $29 < 1.4. As in the previous cases, at each 
deformation the potential energy rises steeply with in- 
creasing A2, as shown in the lower panel of Fig. 7. 

The dynamic path computed in the restricted 2D col- 


FIG. 7. (Color online) Effective collective potential Veg of 
250Fm in the (G20, 830) plane for A2 = 0 (a) and the (20, Az) 
plane for 839 = 0 (b). In each panel energies are normalized 
with respect to the corresponding value of the equilibrium 
minimum, and contours join points on the surface with the 
same energy (in MeV). The energy surfaces are calculated 
with the density functionals DD-PC1 [30], and the pairing 
interaction Eq. (1). 


lective space (329, 330) (dash-dot red curve), and the pro- 
jections of the path determined in the 3D collective space 
(820, 830, A2) (solid curves), are shown in Fig. 8. In the 
(820, 830) plane the 2D and 3D paths are almost indis- 
tinguishable. From the mean-field equilibrium state to 
the fission isomer the reflection-asymmetric shape degree 
of freedom does not contribute, whereas it plays a cru- 
cial role along the dynamic path connecting the isomeric 
state and the outer turning point. The effect of pairing 
fluctuations on the asymmetric fission of ?°°Fm is illus- 
trated in the lower panel of Fig. 8. The dynamic pairing 
correlations are markedly amplified along the fission path 
when the nucleus traverses the inner and outer barriers, 
and this has a significant effect on the effective inertia. 
In fact, even though the 2D and 3D paths are almost 
identical in the plane of axially symmetric shape degrees 
of freedom (620, 830), because of pairing fluctuations in 
the 3D case the fission action integral is reduced by as 
much as ~ 10 units (see Table I). 

The calculated action integrals and resulting fission 
half-lives for 2°4Fm and 7°°Fm are listed in Tab. I. 2D 
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FIG. 8. (Color online) Projections of the 3D dynamic fission 
path (solid curves) of ?°°Fm on the (620, 630) plane for A2 = 
0 (a), and the (620, A2) plane for 839 = 0 (b). The dash- 
dot (red) curve denotes the 2D path computed without the 
inclusion of dynamic pairing correlations. 


TABLE I. Values for the action integral and SF half-lives 
of 7°4Fm and ?°°Fm that correspond to the fission paths 
displayed in Figs. 3 and 8. 


Nucleus Path S(L) logio(Zi/2/yr) 
264m 2D 19.58 —11.03 
3D 14.15 —15.75 
250m 2D 32.09 —0.16 
3D 22.33 —8.64 


denotes the values obtained in two-dimensional collective 
spaces without taking dynamic pairing fluctuations into 
account, whereas 3D labels the values calculated in three- 
dimensional spaces that include the pairing degree of free- 
dom as collective coordinate. For the symmetric sponta- 
neous fission of 7°4Fm the dynamic paths are determined 
in the (820, 822) and (820, B22, A2) collective spaces. The 
path traverses a single fission barrier and the axial sym- 
metry of the fissioning system is fully restored by the 
inclusion of dynamic pairing correlations. As a result, in 
the 3D space the action integral is reduced by ~ 5 units, 


and the predicted half-life is almost five orders of magni- 
tude shorter than in the 2D case that neglects dynamic 
pairing fluctuations. For the asymmetric fissioning nu- 
cleus 7°°Fm we have also shown that, although triaxial 
effects are important in the static case or in the dynamic 
case without pairing fluctuations, the triaxial shape de- 
gree of freedom does not play a role when pairing is in- 
cluded as collective variable, just as in the case of 764Fim. 
Therefore, the SF action integral and fission half-life can 
be determined in the 3D collective space (620, 830, A2). 
Two barriers are traversed by the fission path and, as 
the result of significantly enhanced dynamic pairing cor- 
relations, the action integral is reduced by ~ 10 units as 
compared to the value computed along the 2D path in 
the (820, 830) collective space. The corresponding half- 
life is almost nine orders of magnitude shorter than the 
value predicted in the 2D space (620, (30). 


IV. SUMMARY 


The dynamics of spontaneous fission of 7°4Fm and 
250Fm have been investigated in a theoretical framework 
based on relativistic energy density functionals and, in 
addition to shape deformation degrees of freedom, pair- 
ing correlations have been included as collective coordi- 
nates. Effective potentials and nonperturbative ATD- 
HFB cranking collective inertia tensors have been cal- 
culated using the multidimensionally-constrained rela- 
tivistic mean-field (MDC-RMF) model based on the en- 
ergy density functional DD-PC1, and pairing correlations 
taken into account in the BCS approximation with a sep- 
arable pairing force of finite range. The effect of cou- 
pling between shape and pairing degrees of freedom on 
dynamic (least-action) fission paths, as well as the corre- 
sponding SF half-lives has been analyzed. 

264Fm undergoes symmetric fission into two 132Sn nu- 
clei. Hence, this process can be described in the 3D 
collective space (329, 322, 2), where Az is the Lagrange 
multiplier related to pairing fluctuations via the particle- 
number dispersion operator. The dynamic path that con- 
nects the mean-field ground state and the isomeric state 
of 2°°Fm is also studied in this 3D collective space. For 
both nuclei triaxial deformations reduce the height of 
the static inner barrier by 2 — 3 MeV. However, axial 
symmetry of the fissioning system is fully restored along 
the dynamic paths when pairing is included as collective 
coordinate and, simultaneously, pairing correlations are 
significantly enhanced. The description of asymmetric 
spontaneous fission of ?°°Fm necessitates the inclusion 
of the octupole (reflection-asymmetric) degree of freedom 
639 and, in principle, calculations should be carried out 
in the full 4D space spanned by the collective coordi- 
nates (320, 322, 830, 2). However, since the triaxial de- 
gree of freedom does not play a role in the dynamic case 
that includes pairing fluctuations, the fission action inte- 
gral could be computed along the dynamic path in the 
symmetry-restricted 3D collective space (20, 830, A2). 
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The octupole deformation degree of freedom becomes 
crucial beyond the isomeric state, and pairing correla- 
tions display a pronounced increase when the path tra- 
verses the inner and outer barriers. Consistent with the 
findings of Ref. [21], it has been shown that the inclusion 
of pairing correlations in the space of collective coordi- 
nates, that is, the dynamical coupling between shape and 
pairing degrees of freedom, reduces the fission action inte- 
gral by several units (more than five in the case of 7°4Fm, 
and almost ten for 7°°Fm) and, therefore, has a dramatic 
effect on the calculated SF half-lives. 
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